Interactions and superconductivity in heavily doped M0S2 
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We analyze the microscopic origin and the physical propertis of the superconducting phase recently 
observed in M0S2 . We show how the combination of the valley structure of the conduction band, 
the density dependence of the screening of the long range Coulomb interactions, the short range 
electronic repulsion, and the relative weakness of the electron-phonon interactions, makes possible 
the existence of a phase where the superconducting order parameter has opposite signs in different 
valleys, resembling the superconductivity found in the pnictides and cuprates. 
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Molybdenum disulfide (M0S2) is a layered semicon- 
ductor which can be exfoliated down to monolayer unit 
cells pQ, like graphene [5H3]. The existence of an en- 
ergy gap makes M0S2 a convenient material for nano- 
electronics OE]. Metallic behavior can be induced, also 
like in graphene, by means of electric field effects or by 
doping, and the corresponding Fermi surface is typically 
made up by inequivalent Fermi pockets [7Hll], defining 
a valley degree of freedom which is strongly entangled 
with the spin degree of freedom [T5j, and it can be fur- 
ther controlled and manipulated, opening promising per- 
spectives for spintronics. At high carrier concentrations 
(n ~ 10 14 cm~ 2 ), and in the presence of high-K dielectrics, 
M0S2 has also be shown to undergo a superconducting 
transition, with a doping-dependent critical temperature 
T c (n) which exhibits a maximum as function of n and 
drops to zero at sufficiently large values of n [Ml [14] . 

A ferromagnetic behavior has been also reported in 
MoS 2 [TfjrtTB] , and it has been related to edges or to the 
existence of defects [TS1 [2U]. The magnetic properties 
of M0S2 nanoribbons indicate that the electron-electron 
interactions are non negligible. The combination of sig- 
nificant electron-electron interactions and a two dimen- 
sional Fermi surface made up of many pockets is also 
a hallmark of the cuprate and pnictide superconductors 
|21j . where the superconducting gap has a c£-wave sym- 
metry (cuprates) or opposite sign in different pieces of 
the Fermi surfaces (pnictides) [21]. A related gap struc- 
ture has been also proposed for heavily doped graphene 
when the electron-electron interaction is sufficiently large 

In the present work we study the origin of supercon- 
ductivity in heavily doped M0S2 , by considering the role 
of both electron-electron and electron-phonon interac- 
tions. We analyze first the general features of the ef- 
fective interaction between charge carriers, and we make 
semi-quantitative estimates of the strength of the differ- 
ent contributions to the effective coupling. We discuss 
next the competition between the electron-electron and 
electron-phonon interactions, and the possible types of 
superconductivity that emerge. Although a quantitative 
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FIG. 1: Sketch of the intravalley (a) and intervalley (b) 
scattering processes between Cooper pairs in M0S2. Pan- 
els (c) and (d): corresponding Feynman diagrams associated 
to intra- and intervalley scattering. Full arrows are electron 
propagators, and wavy lines are effective interactions. The 
possible unconventional superconducting phases discussed in 
the text are sketched in panels (e) and (f). 



determination of the superconducting T c is outside the 
scope of our work, the present analysis suggests that su- 
perconductivity in M0S2 is likely to be induced by the 
electron-electron interaction, and that a superconduct- 
ing phase with a non trivial gap structure is possible. 

Effective interactions. Following the experimental re- 
sults of Ref. [131 E] , we assume that the carriers leading 
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to the superconducting phase are electron like, confined 
in the first M0S2 layer closer to the high-dielectric gate. 
The validity of this approximation will be discussed later. 
We also assume that, similar to monolayer M0S2, the 
electron-like carriers are located in the two incquivalent 
Fermi pockets centered at the K and K' corners of the 
Brillouin Zone (see Fig. [I]), having thus a sizable d-orbital 
character with main d^ z 2_ r 2 orbital component. At suf- 
ficiently large concentrations six additional inequivalent 
Q valleys start to be filled, located half way between the 
r and K point, with primary d x 2_ y 2 and d xy Mo orbital 
character. These secondary valleys at higher charge den- 
sity do not play a relevant role in our main discussion 
and will be therefore neglected. Their possible effect will 
be however addressed in the final discussion. 

We explore here the possibility that superconductiv- 
ity is induced by effective electron-electron interactions, 
associated with the direct Coulomb interaction between 
charge carriers or with the effective coupling induced by 
phonons. As in Ref. [22) , we consider singlet super- 
conductivity. Because of the time invariance symmetry, 
the Cooper pair electrons (k f, — k I) reside in differ- 
ent valleys. We can now classify, as sketched in Fig. [T] 
the interactions leading to scattering of the Cooper pairs 
into intra- and intervalley couplings, namely Vi n t ra (q, uS), 
Tenter (q, w), where q and u> are the exchanged momen- 
tum and frequency. We consider only scattering pro- 
cesses of carriers near the Fermi surfaces, which are as- 
sumed to be isotropic and centered at K and K', and we 
neglect the frequency dependence of the interaction. 

The classification of the interaction in an intravalley 
and an intervalley component allows us to define the di- 
mensionless coupling constants [22] 
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where a labels the intra- and inter- valley scattering, and 
p(e F ) — m e f f I \2ttH 2 ) is the density of states at the Fermi 
level per valley and per spin in terms of the effective mass 
m e ff. Typical values are m e ff ~ 0.5mo [231 [21], where 
mo is the free electron mass. For realistic charge con- 
centrations the Fermi wavevector is much smaller than 
the dimensions of the Brillouin Zone, kp <C |K|, so that 
we can in good approximation neglect the momentum 
dependence of V^„ tra (q), T^„ 4e ,.(q), except for electron- 
electron intravalley scattering, as discussed below. 

The existence of superconductivity requires Xmtra ± 
Winter < 0, where the choice of the sign depends on 
whether the gaps in the two valleys have equal or opposite 
signs. The electron-phonon coupling leads to attractive 
interactions, \ e ~ ph < 0, while the electron-electron cou- 
plings lead to repulsive interactions, A e_e > 0. However, 
repulsive interactions can also lead to superconductivity 
provided that A inter > A intra > 0. 

Electron-phonon interaction. The electron-phonon in- 
teraction in single layer M0S2 has been evaluated in Ref. 



[25] . finding relevant contributions from three acoustic 
modes and from six optical modes. Note that the cou- 
pling to the acoustic modes vanishes when the phonon 
wavevector approaches zero. The leading couplings iden- 
tified in Ref. [22] are thus the ones to the polar LO 
modes, Hujlo ~ 0.048 eV, and to the homopolar mode, 
huho ~ 0.05 eV, where the S atoms oscillate out of the 
plane [25] . The homopolar mode contributes to intraval- 
ley scattering, while the LO mode contributes to intraval- 
ley scattering, through the induced electric polarization, 
and also to intervalley scattering. 

Using the notation in [25] , we define an effective in- 
teraction, at frequencies much smaller than the phonon 
energies, as 
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where g^o ~ 0.098 eV is the long wavelength polar cou- 
pling, D LO w 2.6 eV A" 1 and D ho rj 4.1 eV A" 1 are de- 
formation potentials for the LO and homopolar modes, 
respectively, and O is the area of the unit cell. We take 
for M the mass of the sulfur atom, which is much lighter 
than Mo, and thus expected to be dominant. Adding the 
three contributions, we find thus: 



Kntra ~ "0-36, 

\ ph ~ — n i ^ 

A inter ~ 
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which account for the respective intravalley and interval- 
ley electron-phonon coupling constants. 

Electron- electron interaction. After having estimated 
the electron-phonon coupling, we address now the 
electron-electron repulsive interaction. Intravalley scat- 
tering is operative only for small momenta/large dis- 
tance, where, as in Ref. [22]; the electron- electron inter- 
action is determined by the screened Coulomb potential: 
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where eo is the dielectric constant of the environment, 
and qpT = 27re 2 p(e^)/eo is the Thomas-Fermi wave- 
vector. 

On the other hand, the contribution of the electron- 
electron interaction to intervalley scattering is associated 
with the short range part of the Coulomb potential. The 
leading term in this interaction is thus given by the Hub- 
bard term, namely the repulsion between two electrons 
with opposite spin in the same atomic orbital. Since, 
as we discussed above, electronic states close to the K 
and K' points have a dominant Mo 4d character, we can 
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FIG. 2: Intravalley and intervalley dimensionless couplings 
Kntra> Knter due to electron-electron interactions as functions 
of density n, for different values of eo and Uid- Full lines: 
Ai„tra- (from top to bottom: eo = 10,20,30,40,50); broken 
lines: \i„t er (from top to bottom: U44 = 10 eV, 8 eV, 6 eV, 
4 eV, 2 eV). 
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FIG. 3: Relation between the critical carrier density, n c , 
required for the existence of superconductivity and the di- 
electric constant of the environment, eo- Solid blue curve: 
U^d = 2 eV, dashed red curve: Uid = 4 eV. 



therefore approximate: 

VLtr * U id x Q, (5) 

resulting thus in a k-independent and density- 
independent interaction. 

Results. We discuss now the consequences of the 
electron-phonon vs. electron-electron and the intra- vs. 
intervalley interaction in regards to the superconducting 
order. Along this line, unfortunately, the values of the 
dielectric constant, eo, and of the Hubbard term, Uid, are 
not well known. The experiments reported in [131114] are 
done in the presence of a dielectric with a high value of 
eo- The ionization energies of the Mo atom allow us 
to make an order of magnitude estimate of Uid- Alterna- 
tively, we can calculate the Coulomb integral for the Ad 
orbitals of Mo, or replace the orbital by a charged sphere 
with the same radius. All these different approaches con- 
verge on the order of magnitude of the Hubbard repul- 



sion, but a more accurate determination of its specific 
value is lacking. In this situation, we will consider Uid 
and eo as a free variable parameters and we will inves- 
tigate their effects on the superconducting properties of 
MoS 2 in the range e 10 50 and U id m 2 -f- 10 eV. 

Fig. [5J shows the dependence of the intravalley and in- 
tervalley electron-electron coupling constants Xi~ t e ra and 
Knter 0X1 ^ ne specific values of eo and Uid- In particu- 
lar, the intravalley coupling \Xntra shows an initial sig- 
nificant dependence on density, due to the screening in- 
crease with density for qpp < kp, wheeras Xl^tra satu- 
rates for qpp ^> kp. A comparison with the correspond- 
ing coupling constants for electron-phonon interaction, 
provided by Eq. ^ , shows that the values of X\~i ra and 
^inter m the limit of high dielectric constant and mod- 
erate Uid are larger in this regime than those of \Xntra 
and Antral suggesting thus that superconductivity is due 
to the electron-electron interaction. Furthermore, since 
Knter = Knter + Knter > °> tne superconducting phase 
is expected to have gaps with opposite signs in the two 
valleys, as sketched in Fig. [l^. Note that superconduc- 
tivity of this type is possible only when the carrier den- 
sity n satisfies the condition Aj n t ra (n) — \{ n ter < 0. This 
inequality defines a critical density n c above which su- 
perconductivity appears. The dependence of n c on the 
dielectric constant of the environment is shown in Fig. [3] 

Apart from the more likely symmetry of the gap dis- 
cussed above, it is worth to notice that an additional su- 
perconducting phase is possible due to the modulation of 
the electron-electron interaction due to screening. This 
mechanism always leads to an order parameter with a 
modulated fc-dependent gap within each valley [57] (see 
Fig. [jj). However, the resulting critical temperatures 
for this phase are typically very low [25], so that it is 
very unlike to be related to the experimental evidence of 
superconductivity in this material. 

Superconducting phase. The change in sign of the su- 
perconducting gap in the two valleys implies that the 
superconducting phase in M0S2 have unusual proper- 
ties with respect to conventional superconductors. More 
specifically we can remind [57] : i) Elastic scattering is 
pair breaking, leading to the suppression of supercon- 
ductivity when vf/£> |A|, where Dp = Hkp/m e ff is the 
Fermi velocity, I is the elastic mean free path, and A is 
the superconducting gap; ii) Strong scatterers induce lo- 
calized Andreev states within the gap; iii) Andreev states 
also appear at certain edges. However, unlike graphene, 
it should be reminded taht M0S2 also shows a significant 
spin-orbit coupling, so that the combination of a non triv- 
ial superconducting order parameter with the spin-orbit 
coupling can lead to other interesting properties. 

Open questions. An intriguing role, in the above sce- 
nario, is played by the specific value of the effective mass 
m e f f . Apart from the determination of the effective cou- 
plings, such parameter is relevant in assessing the ro- 
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bustness of the Fermi surface structure depicted in Fig. 
[TJ The previous analysis was based on the assumption 
of an electronic structure similar to thant in monolayer 
M0S2, where the carriers occupy two inequivalent valleys 
corresponding at the absolute minima in the conduction 
band at the K and K' points. However, it should be re- 
minded that six inequivalent secondary minima are also 
predicted at the Q points, midway between the Y and the 
K and K' points, shown by the dots in Fig. [T] In bilayer 
and multilayer M0S2 the minima at the Q points are ex- 
pected to lie below the valleys at the K and K' points, 
and in monolayer M0S2 the minima in Q start to be filled 
at sufficiently high carrier densities, once the pockets in 
K and K' are partially filled. Since the relative filling of 
the pockets at K-K' and Q points depends on the cor- 
responding effective masses, an accurate determination 
of m e f f could also assess the possible presence of Fermi 
pockets at the Q points. It should be also noted that the 
value of the effective mass m e f / determines the screening 
properties and, as a result, the distribution of the gate- 
induced charge in the multilayer M0S2. In Ref. [2] it 
was suggested that, because of the strong electric field, 
carriers are mainly localized in the first layer close to the 
high-dielectric gate (note that the high carrier densities 
enhance the screening and thus the confinement). 

Theoretical first-principles-based calculations indicate 
a mass m e ff of the order m e ff ~ 0.4 — lmrj [TJ- 
[TTJ [23J [2H [2H] , whereas early experimental measurements 
suggest m e ff ranging from ~ 0.01to [2t9 to ~ Im [50T 
[32] , In the absence of a definitive estimate of m e ff, we 
checked the dependence of the induced charge distribu- 
tion on m e ff by performing a Thomas- Fermi-like calcu- 
lation using a similar model for the screening of a M0S2 
multilayer as in Ref. [S3]- The results for n « 10 14 cm -2 
using different possible values of m e ft, are shown in 
Fig. [4] They confirm that for m e // > O.lmo the charge 
is indeed confined only in the layer closest to the high-K 
gate [13] . Note that, since the charge is strongly concen- 
trated in the first layer, an appropriate generalization of 
the model on a discrete layer distribution was here em- 
ployed, whereas the continuum original model of Ref. [33] 
would be inappropriate. Our results corroborate thus the 
initial assumption that the single layer model is valid and 
only the K and K' valleys are occupied for the densities 
that lead to superconductivity in [T3J E] . 

It is worth to mention however that, exotic supercon- 
ducting phases, with a gap with different signs in different 
valleys as sketched in Fig. [l^, are expected to appear also 
in the cases when only the Fermi pockets at the Q points 
are filled (expected in the low density regime when the 
carriers are spread among a few layers) or even in the 
case at higher fillings where both pockets (Q and K) are 
occupied. Such exotic phases would share a similar phe- 
nomenology as in the case here considered of only two 
valleys at the K and K 1 points, in particular in regards 
to the features (i)-(iii) mentioned above. 
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FIG. 4: Distribution of the electric charge per layer n[i) in 
multilayer M0S2 under gated conditions for n — I0 14 cm -2 , 
for different values of the effective mass. 



As final consideration, it should be remarked that our 
analysis implies that that the main dependence of the 
superconducting phase on the charge density is by the 
suppression of the long range Coulomb repulsion, due to 
the higher screening at higher charge concentrations (see 
Fig. [2| . Hence, our analysis at this level does not account 
for the decrease of T c observed at high carrier densities 
n > 1.3 x 10 14 cm~ 2 [14]. A possible explanation for this 
behavior could be the change in the orbital character 
of the states close to K and K' point, which loose Ad 
orbital character at higher doping, resulting in a weaker 
Hubbard-like intervalley interaction, or the change in the 
topology of the Fermi surface, as the Q valleys start to 
be filled. 

Conclusions. In this paper we have analyzed the ap- 
pearing of a superconducting phase of M0S2 at high car- 
rier concentrations and for strong screening of the long 
range Coulomb potential. The significant short range re- 
pulsion between carriers at the conduction band allows 
for a superconducting phase induced by the electron- 
electron interaction, where the gap acquires opposite 
signs in the two inequivalent pockets of the conduction 
band. This superconducting state, similar to that found 
in the cuprate and pnictide superconductors, is expected 
to show interesting topological features, such as Andreev 
states at edges and grain boundaries. 
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